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This  stud/  undertakes  to  determine  the  feasibility  of  using  present  do/  stropdown  inertial 
navigation  systems  to  motion  compensate  synthetic  aperture  radars,  and  In  particular,  motion 
compensation  using  the  autofocus  technique.  It  is  important  to  understand  that  the  autofocus 
technique  is  only  of  many  methods  being  investigated  in  the  field  as  the  best  means  to  motion 
compensate  a  highly  accurate  rador  detection  system.  It  would  be  beyond  the  capabilities  of  myself 
and  that  of  a  thesis  to  undertake  a  thorough  study  of  all  the  motton  compensation  techniques 
available.  I  have  attempted  to  make  the  reader  understand  this  particular  form  of  motion 
compensation  and  its  very  reel  possibilitas  as  the  most  capable  for  highly  dynamic  aircraft.  The 
ability  to  use  synthetic  aperture  radar  in  extremely  maneuverable  aircraft  is  only  now  becoming 
possible  due  to  very  small  but  highly  capable  onboard  digital  computers.  The  combination  of  these 
computers  and  the  robustness  of  the  strapdown  inertial  systems  and  synthetic  aperture  radars 
promise  to  make  the  aircraft  of  tomorrow  a  formidable  weapon  system. 

It  was  very  interesting  to  me  in  my  capacity  as  an  F-lil  Weapon  Systems  Officer  to 
investigate  the  possibilities  Inherent  in  synthetic  aperture  radar  and  strapdown  inertial 
navigation  systems.  Their  low  life-cycle-cost  and  maintainablity  seem  to  make  it  impcritive  that 
more  research  is  undertaken  in  this  field. 

I  would  like  to  take  this  opportunity  to  thank  those  involved  in  my  quest  for  knowledge  at 
AFIT.  My  thesis  advisor,  Dr.  Slpurls,  was  extremely  helpful  and  always  ready  to  steer  me  in  the 
right  direction  if  I  happened  to  stray  from  the  straight  and  narrow.  Lt  Col  Watt,  my  faculty 
advisor ,  was  very  helpful  when  I  thought  the  going  was  getting  too  tough  and  encouraged  me  to  keep 
at  it  just  a  little  longer.  And  finally,  my  wife  and  children,  who  always  kept  things  in  perspective 
and  made  the  experience  worthwhile. 
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AUTOFOCUS  MOTION  COMPENSATION  FOR  SYNTHETIC  APERTURE  RADAR 
AND  ITS  COMPATIBILITY  WITH  STRAPDOWN  iNERTIAL  NAVIGATION 
SENSORS  ON  HIGHLY  MANEUVERABLE  AIRCRAFT 

I.  Introduction 

Background 

As  modern  military  aircraft  hsve  become  more  sophisticated,  aid  more  expensive,  their 
ability  to  pinpoint  targets  on  the  ground  and  strike  those  targets  with  weapons  equally 
sophisticated  and  expensive  has  lagged  behind  In  the  past,  gimbaled  Inertial  Navigation  Systems 
(INS)  have  provided  precise  autonomous  navigation  and  position  information  for  military 
aircraft.  The  mechanical  complexity  and  penalties  of  gimbaled  systems  such  as  maintainability, 
reliability,  and  Life-Cycle-Cost  (LCC)  limit  their  application.  These  systems  cannot  realize  the 
low  cost  and  reliablity  of  strapdown  inertial  systems  utilizing  ring  laser  gyroscopes  or  tuned 
gyroscopes  such  as  the  SPN-GEANS  or  N-73.  Strapdown  Inertial  Navigation  Systems  utilizing 

l 

either  tuned  gyroscopes  or  ring  laser  gyroscopes  and  advanced  accelerometers  mey  alleviate  the 
poor  maintenance  record  of  inertial  platforms  which  have  plagued  the  advanced  aircraft  of  today. 
These  strapdown  inertial  systems  must  significantly  reduce  their  LCC  without  sacrificing 
performance  and  accuracy,  it  Is  felt  that  with  the  availability  of  small  high-speed  digital 
minicomputers  the  strapdown  INS  has  finally  become  an  effective  system.  Also,  strapdown 
inertial  navigation  systems  have  become  nearly  as  aerate  as  their  predecessors,  and  as  this 
study  will  show,  capable  af  handling  the  motion  compensation  tasks  of  modern  radar  applications. 

Problem 

Modern  military  aircraft's  ability  to  navigate  using  their  onboard  radar  to  update  the 
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inertial  navigation  system  has  been  held  back  by  the  ,jck  of  high  definition  information  caning 
from  the  radar  system.  Although  the  increasing  sophistication  of  modem  military  aircraft  has 
made  h  possible  for  them  to  handle  much  more  difficult  tasks  and  on  ever  expanding  role  in  the  Air 
Order  of  Battle,  their  ability  to  seek  out  and  identify  targets  more  accurately  has  not  kept  up  with 
other  more  technologically  advanced  areas  of  aviation.  Synthetic  Aperture  Radar  (SAR)  me/  be 
able  to  provide  a  much  bettor  picture  of  the  terrain  with  which  to  update  the  onboard  navigation 
systems.  Synthetic  Aperture  Radar's  promise  in  the  area  of  enhanced  resolution  of  the  radar 
picture  me/  be  held  back  only  by  Its  lack  of  highly  accurate  Information  from  the  Inertial 
navigation  system  with  which  to  allow  motion  compensation.  The  combination  of  improved 
navigation  and  target  tracking  systems  promise  to  make  the  aircraft  of  tomorrow  far  superior  to 
anything  now  flying. 

SCQQfi 

It  is  the  intention  of  this  study  to  show  that  the  motion  compensation  requirements  of 
synthetic  aperture  radars  aboard  these  modern  military  aircraft  can  be  met  by  strapdown  inertial 
navigation  systems  that  are  available  at  present.  Fundament''1  to  motion  compensation  and  high 
resolution  Synthetic  Aperture  Radar  is  the  concept  of  focusing.  In  particular,  the  autofocus 
technique  of  motion  compensation  is  explored.  The  autofccus  technique  is  investigated  as  a  means 
of  expanding  the  range-squint  angle  boundaries  of  synthetic  aperture  radar  operation. 

Development 

The  differences  in  Reel  and  Synthetic  Aperture  Radar  are  first  discussed  giving  the  reader  a 
refresher  in  radar  theory  and  leading  to  discussion  of  aperture  length  and  integration  time 
requirements  as  well  as  SAR  focusing  end  aperture  weighting  SAR  phase  errors  and  image  quality 
will  then  be  covered  with  particular  attention  being  given  to  phase er*urs  caused  by  SAR  aircraft 


motion  measurement  errors  and  tneir  visual  effects.  Unaided  INS  error  characteristics  will  then 


be  explored  showing  the  feasibility  of  a  one  nautical  mile  per  hour  INS  for  performing  the  motion 
compensation  requirements.  Then  oided  INS  performance  capabilities  are  demonstrated  for  typical 
strapdown  Inertial  Navigation  Systems  such  as  the  SPN-0EANS,  N-73,or  Ring  Laser  Oyros.  The 
previour  sections  have  merely  been  designed  to  show  the  conceptual  compatibility  between 
strapdown  INS  performance  and  SAR  motion  compensation  requirements.  Next,  SAR  motion 
compensation  requirements  are  covered  to  illustrate  ttie  INS  performance  necessary  to  achieve  the 
quadratic  phase  error  required  for  SAR  focusing  for  specific  operating  conditions.  Once  the 
compatibility  between  INS  performance  and  SAR  motion  compensation  requirements  have  been 
shown,  it  is  still  necsssary  to  use  some  sort  of  focusing  technique  to  obtain  high  resolution  images 
under  more  severe  aircraft  maneuvering.  This  study  is  concluded  with  an  explanation  of  the 
Autofocus  technique  for  motion  compensation  giving  SAR  performance  envelopes  with  and  without 
the  Autofocus  process  being  used. 


il.  Real  and  Synthetic  Rnrinr 


Real  array  imaging  radars  obtain  y-direction  resolution  by  virtue  of  the  real  antenna 
beam  coverage  on  the  ground  which  is  determined  by  the  physical  size  of  the  antenna  The 
y-direction  being  the  along  trade  or  direction  in  which  the  aircraft  is  moving.  Finer  resolution 
will  require  longer  antennas,  which  because  of  their  si2e  can  only  be  .nounted  lengthwise  along 
the  aircraft.  This  results  in  an  antenna  beam  looking  to  the  side  of  the  aircraft.  As  will  be  seen 
later,  incremental  Doppler  shift  of  adjacent  resolution  cells  in  the  y-direction  can  only  be 
seperated  if  the  antenna  is  looking  to  one  side  rather  than  straight  ahead.  This  is  typically  the 
method  used  in  reconnaisance  aircraft  and  not  useful  in  an  aircraft  engaged  in  weapon  delivery 
or  navigation  updating  with  which  this  study  is  concerned.  To  make  synthetic  radars  useful  in 
tactical  applications,  the  radar  will  have  to  look  slightly  to  the  side  rather  than  straight  ahead. 
This  angle  between  the  direction  the  aircraft  is  moving  and  the  direction  in  which  the  radar  is 
looking  will  be  called  the  squint  angle. 

Figure  2. 1  shows  a  simplified  diagram  of  the  geometry  of  a  real  array  imaging  radar.  The 
aircraft  flies  along  the  y-direction  at  velocity  V.  Its  long  antenna  produces  a  fan  beam  which 
illuminates  the  ground  below.  For  simplicities  sake,  the  squint  angle  is  taken  to  be  ninety 
degrees.  It  is  Important  to  note  that  the  beams  width  is  narrower  in  the  y-direction  than  its 
length  in  the  x-direction.  The  y-direetjon  resolution  is  determined  by  the  beamwidth  while  the 
x-direction  resolution  is  determined  by  the  pulse  length  and  is  independent  of  the  beamwidth.  In 
an  aircraft  using  real  imaging  radar ,  the  antenna  alternately  transmits  end  receives  energy  and 
displays  the  ground  returns  received  on  a  cathodB  ray  tube  or  records  the  images  on  film.  The 
resolution  needed  to  determine  objects  will  be  covered  first 


Fig.  2.!  Geometry  of  a  Real  Array  Imaging  Radar  (17:33) 


Rfisniirtinn  Call  Size 

It  3tou)d  be  helpful  to  discuss  required  resolution  for  determining  objects  before 
covering  achievable  radar  mapping  resolution.  Table  2. 1  gives  the  required  resolution  cell 
sizes  for  various  objects.  Large  objects,  such  as  coastlines,  cities,  and  mountains,  can  be 
recognized  with  150m  resolution  elements.  Vehicles,  houses,  and  buildings  however  need 
resolution  elements  along  the  order  of  3  meters.  This  gives  a  reference  for  determining  the 
resolution  required  to  go  after  certain  targets  using  a  tactical  aircraft 


TABLE  2.1  RESOLUTION  CELL  SIZE  FOR  VARIOUS  OBJECTS  ( 1 3:3) 


ITEM 

SQUARE  CELL  SIZE  (meters) 

Coastlines,  Cities,  Mountains 

150 

Major  Highways,  Large  Airfields 

30 

City  Streets,  Large  Buildings 

15 

Vehicles,  Houses,  Buildings 

3, 

X-Directinn  Resolution 

Pulse  radars  measure  range  by  transmitting  a  pulse  and  timing  the  returned  pulse  from 


Fig.  2.3  Resolution  in  the  X- Direction  (13:5) 


F3 


If  the  range  gates  are  of  duration  T ,  two  adjacent  points  on  the  ground,  A  and  B ,  which  are 
aRx  apart,  can  be  resolved  by  their  returns  occuring  in  two  adjacent  range  gates  as  shown. 

Electromagnetic  energy  travels  at  the  speed  of  light,  c,  so  the  time  for  the  energy  to  travel  to  the 
target  and  return  is  given  by 


2R  .  t  +T  _  2(R  *  AR) 
c  '  c 


(sec)  (20 


Subtracting  equations  (2.1 )  gives 


T  = 


2AR 


(sec)  (2.2) 


and  from  observation  of  Figure  2.3 


ARy  s 


AR 


cos  a  2cos  a 


(m)  (2.3) 


where  aRx  is  the  resolution  in  the  x-direction  or  perpendicular  to  the  ground  track.  It  can  be 
seen  from  equation  (2.3)  that  resolution  in  the  x-direction  is  directly  proportional  to  the 
pulsewidth  X.  Therefore,  shorter  pulsewidths  produce,  better  resolution.  In  radars  which  are 
matched,  i.e.  the  product  of  pulse  duration  X  and  the  receiver  bamfcridth  W  is  unity,  one  gets 

cx  cl tv)  c  ,  .  , 

ARX  - - =  — =  -  (m)  (2.4) 

2cos  a  2cos  <x  2Wcos  a 


In  the  pulse  radar  shown  in  Figure  2.2,  the  maximum  measurable  range  R^  is  a  function  of 
interpulse  period  T ,  and  is  given  by 
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(2.5) 


^max  =  ~  (Km) 

If  R  is  beyond  this  .‘ange,  the  elapsed  time  between  transmitted  pulse  and  target  return  pulse 
will  include  multiples  of  intsrpulse  period  T ,  making  the  measured  range  ambiguous. 

Y-Direction  Resolution 

In  real  arrav  Imaging  radars  resolution  in  the  y-direction  is  proportional  to  the  antenna 
beamwidth.  A  rectangular  antenna  with  uniform  current  distribution  will  produce  a  (sin  x)/x 
antenna  pattern  (4:385-389,  12:2-20).  The  approximate  expressions  for  half-power 
beamwidth  of  this  antenna  are 

T|  =--A  c=52A  (decrees)  (2  6) 

'  '2 

where  Tl  and  e  are  half-power  beamwidths  in  degrees,  A  is  the  wavelength  and  1  ,and  12  are  the 

sides  of  the  antenna  A  circular  antenna  with  uniform  current  distribution  will  produce  a 
radiation  field  intensity  involving  first  order  Bessel  functions  ( 1 2:2-24).  The  half-power 
beamwidth  in  this  case  will  be 

y  =  (deqrees)  »  —  (radians)  (2.7) 

0  O' 

where  y  is  the  half- power  beamwidth  and  D  is  the  diameter  of  the  antenna  In  further 

discussion  the  value  of  y  will  be  approximated  by  y*A/D  where  y  is  given  in  radians.  It  is 
now  easy  to  determine  the  resolution  in  the  y-dtrection.  From  Figure  2. 1  the  resolution  in  the 
y-direction  is 
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r.  m 


: 


ARg  =  Rti  (m)  (2.8) 

where  ARy  is  the  resolution  in  the  y-direction,  R  is  the  range  to  the  resolution  cell  and  T|  is 

ahtehna  azimuth  beamwidth.  If  the  value  of  Tt  is  approximated  by  X/ 1 ,  .the  resolution  in  the 
y-direction  is  given  by 

aRu  =i-A  (m)  (2.9) 

As  seen  by  equation  ( 2.9),  the  resolution  in  the  y-direction  is  dependent  on  range  R.  Therefore 
as  the  range  gets  smaller,  resolution  improves.  This  property  in  reel  array  imaging  radars 
make  ground  points  closer  to  the  aircraft  resolvable  while  points  further  awey  merge  together.  , 
In  both  real  and  synthetic  array  imaging  radars,  resolution  in  the  x-direction  depends  on  the 
time  of  arrival  of  the  target  return  signal.  Because  of  this  phenomenon,  equidistance  points 
from  the  radar  will  all  converge  to  the  same  range  gste  and  these  points  will  all  occur  at  the 
same  x-resolution  cell.  To  someone  at  the  position  of  the  aircraft,  these  points  will  show  up  at 
three  distinct  locations. 

The  need  for  ground  maps  without  using  large  antennas  led  to  synthetic  array  radars. 
These  radars  use  pulse  duration,  as  in  the  case  of  reel  array  imaging  radars,  for  x-direction 
resolution.  However ,  y-direction  resolution  is  determined  by  using  Doppler  shifts  or  spectral 
analysis  as  opposed  to  beamwidth.  Therefore,  the  long  antennna  producing  a  narrow  beam  is 
replaced  in  the  synthetic  radar  by  a  great  deal  of  spectral  signal  processing.  By  using  synthetic 
signal  processing  techniques  one  can  obtain  high  resolution  ground  maps  using  standard  antennas 
in  aircraft.  Since  the  basis  of  resolution  in  the  y-direction  is  frequency  shift,  which  is  the  rate 
of  change  of  phase,  a  discussion  of  these  topics  will  be  given. 
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If  u  iransmitted  waveform  given  by 


Vt  =  sir.  2irf0t  (m/sec)'  (2.10) 

Is  emitted  from  an  antenna  and  this  sijyial  is  reflected  from  a  moving  target  with  radar-target 
range  R  given  by  the  equation 

R  =  R0  +  Rt  (km)  (2.11) 

where  R0  is  the  initial  range,  R  Is  target's  closing  r8te  in  the  direction  of  R,  and  t  is  the  time. 

The  returned  signal  at  time  t  will  be  the  signal  which  was  transmitted  &t  seconds  ago  with,  At 
given  by  t  in  equation  (2.1).  If  one  substitutes  t=t-*t  in  ( 2. 1 0)  the  received  signal  is 


(m/sec) 


(2.12) 


Once  terms  involving  t  are  grouped  together 


♦  -  2(2irf0)Ro 
c 


(2.13) 


Notice  that  the  coefficient  of  time  t  in  the  argument  represents  frequency  while  the  part  not 
involving  time  represents  a  phase  The  frequency  and  phaee  of  the  returned  signal  will  be 


2tT fr  =  2TT f 0  -  I1 ?.(2_ ,°)  ,  a.  =  -  2(2Trf°>o  (2. 14) 

f'  ■  A 


To  obtain  the  difference  of  transmit  and  receive  frequencies,  fP-f0,  equation  (2.14)  is 
rearranged  to  arr  ive  at 

fd  =  fr-fo  =  l±j-  =  lil  (Hz)  (2.15) 

c  c/f©  A 

where  fd  is  recognized  as  the  doppler  shift  and  c/f0  is  the  wave  length  A.  It  can  be  seen  that  the 

return  from  a  moving  target  will  be  shifted  in  frequency  directly  proportional  to  radar-target 
closing  rate  and  inversely  proportional  to  wavelength.  Also,  the  phase  angle  of  the  returned 

signal  $r  as  given  by  (2.14)  can  be  written  as 

*r ='4ir  (rod,ons)  (216) 

Equation  2. 1 6  represents  the  number  of  equivalent  wavelengths  to  the  round  trip  distance  2Rg. 
The  factor  of  2tr  converts  the  number  of  wavelengths  to  angles  in  radians.  Therefore,  the  phass 
change  is  equivalent  to  a  distance;  Similarly,  the  Doppler  shift  of  equation  (2.15)  for  an 
incoming  target  can  be  interpreted  as  the  number  of  wavelengths/ second  that  the  transmitted 
signal  is  pushed  in  due  to  the  closing  rate  R.  The  number  of  wavelengths/second  is  recognized  to 
be  frequency.  To  summarize,  it  should  be  noted  that  phase  is  associated  with  range  while  the 
derivative  of  phase  with  respect  to  time  ( frequency)  is  associated  with  range  rate. 

Again,  resolution  in  the  y-dlrection  Is  looked  at  Using  Figure  2-4,  an  aircraft  flying  at 
velocity  V  parallel  to  the  y-axis  attempts  to  differentiate  between  points  A  and  B.  To  do  this,  the 


It 


frequency  difference  between  the  two  ground  points  A  and  B,  which  are  a  distance  d  apart.  Is 
obtained  Using  the  previous  aquations  and  Figure  2  A  the  Doppler  shifts  for  the  points  are 
extracted 

f  _2Vcos8  f  .  2Vcos(8+&0)  (H7)  (2  17) 

B  _  A  ~ 

If  cos(  0  ♦  A0)  is  expanded  in  terms  of  sines  and  cosines  and  using  the  small  angle  relationships 
of  cosA0  s  1  and  sinA0  s  A0,  the  differences  in  frequency  are 

(Hz)  (2.18) 

.  (2.19) 

Af.=  ^3in0  (Hz).  (2.20) 

1  AR 


Af.=  fQ-  fA=2VA09lf)0 
’  8  A 


From  Figure  2 A  A0  can  be  approximated  by 


R 


where  R  is  the  range.  If  equation  (2. 19)  is  used  in  (2. 18) 


t2 


Therefore  resolution  in  the  y-direction  becomes 


AR  =  d  =  -il  (m)  (2.21) 

y  2Vsin0 

The  resolution  attained  above  is  independent  of  antenna  beemwidth  ( as  was  not  the  case  in 
real  array  imaging  radBrs)  and  depends  primarily  on  attainable  frequency  resolution  *f,  of  the 
signal  processing  equipment  Again  for  a  matched  transmit/receive  system  the  duration  of 
target  illumination  T  should  be  \/iv  What  this  means  is  that  during  the  data  aquisition  time  T 

the  aircraft  flies  a  length  of  L=VT  where  V  is  its  velocity.  Thus,  it  can  be  seen  the  equivalent 
antenna  array  length  becomes 

L  =  VT=^=W^9  <"•>  (2.22) 

Using  the  value  of  L  from  ( 2.22)  another  value  of  y-direction  resolution  is  obtained 

iR*  ■  d  ■  A  (m)  (Z23) 

This  value  of  y-direction  resolution  Is  the  best  attainable  and  is  used  as  a  measure  of  merit  in 
further  calculations.  In  actual  practice,  the  y-direction  resolution  is  obtained  by  a  set  of  tuned 
filters  as  shown  in  Figure  2.5.  This  figure  shows  a  set  of  filters  with  a(sinx)/x  drop-off 
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Fig.  2.5  Actual  and  Theoretical  Y-Direction  Resolution  (13:20) 


characteristics.  The  3-dB  bandwidths  of  the  filters  in  Fig.  2.5  are  approximately  one-half  of 
the  null-to-null  bandwidth.  Referring  to  Fig.  2.5,  a  typical  target  will,  fall  into  two  adjacent 
filters,  eg.,  filters  1  and  2.  Averaging  will  put  the  target  in  the  center  of  these  filters  as 
shown.  To  obtain  proper  resolution  seperatton  filter  number  3  should  be  skipped  and  the  next 
resolvable  target  will  probably  appear  in  filters  4  and  5.  Averaging  will  put  the  target  between 
these  filters  as  shown.  Thus,  it  can  bo  seen  that  the  practical  attainable  resolution  is  3d 

doing  back  to  equation  ( 2.22),  representative  SAR  integration  times  and  aperture  lengths 
can  be  determined  as  shown  in  Table  2.2. 


Table  2.2  Representative  Integration  Times  ami  Aperture  Lengths 
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III.  SAR  Focusing  and  Aperture  Weighting 

The  concept  of  focusing  is  fundamental  to  motion  compensation  and  high  resolution  SAR. 
There  are  several  analytical  approaches  which  are  possible  when  considering  the  topic  of 
focusing.  Among  these  are  vector  addition,  cross  correlation,  and  matched  filtering.  All  these 
techniques  give  approximately  the  same  results.  The  approech  which  will  be  followed  here  is 
that  of  real  beam  antennna  theory  in  which  the  field  intensity  pattern  of  the  antenna  is  related  to 
the  current  distribution  across  the  aperture  by  the  Fourier  transform.  If  the  parabolic 
reflector  real  beam  antenna  of  Figure  3.  f  is  considered,  it  can  be  seen  that  any  ray  from  the 
focus  is  reflected  in  a  direction  parallel  to  the  exis  of  the  parabola  The  distance  traveled  by  any 


Fig.  3.1  Parabolic  Reflector 


ray  from  the  focus  to  the  parabola  and  by  reflection  to  a  plane  perpendicular  to  the  parabola  axis 
is  independent  of  its  path.  This  characteristic  results  in  a  point  source  of  energy  located  at  the 
focus  being  converted  into  a  plane  wave  of  uniform  phaaei  Note  the  similarity  of  Figure  3.2  to 
that  of  Figure  3.1.  A  circular  flight  path  with  its  origin  at  the  map  point  would  also  result  in  a 
SAR  aperture  distribution  having  uniform  phase,  if  it  were  possible  to  maintain  this  flight  path 
the  SAR  aperture  would  be  focused  and  no  motiuon  compensation  would  be  necessary.  This  type 
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Fig.  3.2  Uniform  Phase  Aircraft  Flight  Profile 


of  flight  path  is  impossiole  to  maintain  and  is  in  fact  tactically  unacceptable,  therefore  motion 
compensation  must  be  used  to  compensate  the  phase  of  each  transmitted  radar  pulse  such  that  an 
effective  synthetic  aperture  is  reli2ed  having  uniform  phase  across  its  length.  The  electric  field 

intensity  pattern,  E(  $),  can  be  obtained  as  the  Fourier  transform  of  the  current  distribution 
across  the  aperture.  This  pattern  can  be  determined  for  an  aperture  distribution  of  uniform 
phase  and  constant  amplitude  as 

1/2  L/2 

f  J2rrY  31  n  f  i2Tff  sin  ft 

E(4>)  =  j  A(Y)e  >■  dY  =  A,|  e  ^  dY  (3.1) 

-L/2  -L/2 

...  _  Sint  ff(L/Msin<t>  ]  _  sin  x  ,7 

jr(LA)sinft  "  x 


where  L  is  the  length  of  the  aperture  in  the  Y  direction  and  ft  is  the  azimuth  angle  measured 
from  a  normal  to  the  Y  direction.  If  the  square  of  E(  i})  is  plotted,  the  power  radiation  pattern 


would  be  obtained  as  shown  in  Figure  3.3.  This  should  be  familiar  as  the  (sin  x)/x  pattern  . 
This  pattern  is  applicable  for  both  the  real  beam  parabolic  antenna  of  Figure  3.1  and  the 
synthetic  aperture  antenna  of  F igure  3.2.  Resolution  of  either  antenna  is  inversely  proportional 
to  the  aperture's  length.  It  should  be  noted  that  the  radiation  pattern  corresponding  to 
uniform 


amplitude  distributions  has  a  peak  side-to-main  lobe  ratio  of  -13.2  db.  These  high  side  lobes 
will  adversely  affect  the  quality  of  SAR  imagery  by  filling  in  low  radar  return  terrain,  such  as 
roads,  rivers,  lakes,  and  runways,  with  energy  coming  in  through  the  side  lobes  from 
potentially  high  radar  back  scattering  terrain  such  as  trees,  rocks,  and  buildings.  This 
characteristic  has  the  effect  of  reducing  contrast  or  washing  out  the  SAR  map.  This  problem  can 
be  improved  by  reducing  the  side  lobe  level  through  tapering  or  weighting  the  aperture 
amplitude  distribution  to  decrease  as  a  function  of  distance  from  the  center.  However  the 
reduction  of  the  sidelobes  increases  the  width  of  the  main  loba. 
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IV.  SAR  Phase  Frrors  and  I  may  Quality 


The  ability  to  acheive  a  perfectly  uniform  phase  distribution  across  the  aperture  is 
unrealistic  In  practice  due  to  the  many  ways  in  which  phase  errors  can  occur  in  SAR  systems; 
atmospheric  propagation,  transmitter  nonlinearities,  instabilities,  motion  compensation,  etc. 
The  phase  errors  which  occur  affect  system  performance  by  dagrading  the  antenna  radiation 
pattern,  and  more  generally,  the  system  impulse  response  which  additionally  includes  the  range 
dimension.  As  phase  errors  increase,  the  system  impulse  response  is  degraded  in  many  ways; 
the  siuslobe  levels  are  increased,  the  main  lobe  width  is  increased  resulting  in  a  loss  of 
resolution,  and  the  mainlobe  peak  is  reduced.  This  study  is  primarily  interested  in  phase  errors 
caused  by  aircraft  INS  motion  measurement  errors.  The  effects  produced  by  these  phase  errors 
are  most  easily  analyzed  by  separating  them  into  two  types;  law  and  high  frequency.  Figure  4. 1 
depicts  these  errors. 


The  low  frequency  errors  exhibit  less  than  one  cycle  varietion  over  the  processing 
aperture  and  can  be  represented  as  a  Taylor  series 


(4.1) 


4/(1)  =  4ffR(t) 

A 

xWt )  -  4H  RCo)  +  t  +  jd  r\\  + 

w  (at2/ 2! 

where:  R(0)=the  SAR  aircraft- to- map  point  range  at  the  beginning  of  the  coherent 

integration  time  interval. 

t  a  the  time  referenced  to  the  start  of  coherent  integration. 

The  error  In  the  first  term  of  Eq.  4.2  which  is  time  independent  must  be  small  enough  to  permit 
illumination  of  the  proper  area  to  be  mapped,  cue  the  target,  and  meet  the  weapon  delivery 
requirements.  The  error  in  the  second  term  of  Eq.  4.2  which  is  linear  in  time  and  distance 
results  in  an  azimuth  shift  of  the  system  impulse  response  but  does  not  otherwise  degrade  it.  In 
the  radar  “open- loop"  the  SAR  map  Is  tracked  in  both  azimuth  angle  and  Doppler  and  therefore 
the  magnitude  of  the  second  term  must  be  small  enough  to  insure  that  the  antenna  illuminates  the 
same  ground  patch  that  the  signal  processor  is  filtering.  The  error  in  the  third  term  of  Eq.  4.2 
is  by  far  the  most  important  and  generally  produces  all  the  degrading  effects  described  earlier. 
These  effects  are  shown  In  f igure  4.2  where  the  phase  error  is  from  the  center  to  the  edge  of  the 
aperture.  The  higher  order  low  frequency  phase  errors  affect  the  system  impulse  response  in 
the  same  manner  as  the  quadratic  phase  error  (OPE)  but  their  magnitudes  are  greatly  reduced 
and  do  not  limit  system  performance. 
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Fig.  4.2  Effect  of  Quadratic  Phase  Error  on  Resolution  and  Gain  ( 1 6:9) 


High  Frequency  Errors 

The  effects  of  high  frequency  sinusoidal  phase  errors  on  a  SAR  impulse  function  response 
has  been  known  for  sane  time  under  the  theory  of  paired  echos.  The  following  expression  is  a 

good  approximation  to  the  azimuth  channel  response  v(x)  of  a  SAR  system  for  a  peak  sinusoidal 
phase  error  of  one  radian  or  less.  i 


v(x)  =  J0(*)v0(x)  +  jj,(4/) 


vo(x-7~)  +  i/o(x  + 


f>] 


(4.3) 


where  xp  is  the  phase  error  amplitude,  (J  is  the 


phase  error  radian  frequency,  Jm  is  the  Bessel 


function  of  the  first  kind  and  of  order  m,  V0(x)  is  the  desired  response,  x  is  the  azimuth 


distance,  and  k=4if /XR.  The  desired  response  isj  reduced  in  amplitude  by  the  amount  J0(\jf )  and 
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i$  surrounded  with  a  pair  of  echoes  with  displacements  of  (*)/k  as  shown  in  Figure  4.3. 


Since  phase-corrupted  and  phase-error- free  outputs  have  equal  energies,  energy  is  removed 
from  the  main  lobe  and  redistributed  in  the  sidelobes.  Increase  In  the  side  lobe  levels  is 
undesirable  since  it  leads  to  spurious  target  responses  and  poor  image  contrast  The  peak  side 
lobe  to  main  lobe  ratio  ( PSLR)  is  given  by: 


PSLR  = 


Jjty) 

2 

1 

1 - 

t 

* _ 

Jo(*> 

"  4 

1 

.  2  32  . 

(4.4) 


The  primary  effect  from  high  frequency  random  phase  errors  is  also  shown  in  Figure  4.3  and  is 
to  take  main  lobe  energy  and  distribute  it  over  several  side  lobes  rather  than  concentrating  the 
energy  as  sinusoidal  phase  errors  do.  An  analysis  of  the  integrated  side  lobe  energy  shows  that 
its  ratio  relative  to  the  main  lobe  (ISLR)  is  approximated  by  the  variance  of  high  frequency 
phase  errors  for  errors  a  third  of  a  radian  or  less: 


(4.5) 


Summary 

To  tie  this  section  on  SAR  phase  errors  together,  their  qualitative  effect  on  image  quality 
is  considered  and  can  be  shown  via  Figure  4.4.  Quadratic  phase  error  (QPE)  causes  broadening 
of  the  image  processing  system's  impulse  response  width.  For  90  degress  of  QPE  measured  from 
the  center  of  a  coherent  array  to  its  end  (sometimes  called  the  Rayleigh  limit),  the  main  lobe 


Fig.  4.4  Visual  Effects  of  Poor  QPE,  PSLR,  and  ISLR  ( 1 6:1 2) 


width  is  increased  by  about  8  percent  for  a  uniformly  weighted  array,  this  broadening  causes 
an  effective  decrease  in  image  resolution  and  tends  to  smear  targets.  However,  the  effect  of  QPE 
is  not  limited  to  the  main  lobe.  One  of  its  first  effects  is  to  fill  in  side  lobe  nulls  and  increase 
ISLR  causing  some  of  the  3ame  visual  effects  tnat  ISLR  does.  When  error  levels  are  sufficiently 
high,  PSLR  problems  are  probably  the  easiest  to  observe  and  identify.  Multiple  images  appear 
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in  the  map  as  if  various  shifted  versions  of  the  true  scene  were  overlaid  on  top  of  one  another. 
Increased  PSIR  is  also  accompanied  by  degraded  ISLR  but  (SLR  problems  are  more  often  the 
result  of  phenomena  more  distributed  in  frequency  which  causes  an  overall  increase  in  side  lobe 
levels  without  distinctive  side  lobe  peaks.  This  side  lobe  corruption  may  be  so  severe  that  it 
obliterates  edges  and  contrast  boundaries  in  the  scene.  The  image  appears  washed  out. 


V.  INS  Characteristics  and  Capabilities 


Navigation  by  the  process  of  determining  position  from  self  contained  inertial 
measurements  made  on  board  the  aircraft  is  termed  inertial  navigation.  A  typical  Inertial 
Navigation  System  consists  of  gyros,  accelerometers,  a  computer,  and  associated  equipment.  The 
gyros  are  used  to  establish  an  inertial  reference  frame.  The  components  of  inertial  acceleration 
as  measured  in  this  reference  frame  are  determined  from  the  outputs  of  accelerometers  after 
accounting  for  the  effect  of  gravity  an  the  measurements.  Then  a  single  integration  gives  the 
velocity  and  another  integration  gives  the  position  of  the  aircraft  It  can  be  demonstrated  that 
the  propagation  of  errors  in  an  INS  is  very  accurately  dwribed  by  a  set  of  linear  differential 
equations  under  a  very  broad  set  of  assumptions  (26).  The  minimum  number  of  equations 
necessary  to  simulate  the  INS  errors  is  nine;  three  each  for  velocity,  position,  and  attitude 
errors.  These  are  the  familiar  Pinson  error  states.  If  additional  INS  hardware  and  various 
computational  errors  are  to  be  simulated,  the  number  of  simultaneous  error  equations  required 
increases  on  a  one- for -one  basis. 

Unaided  INS 

The  typical  military  aircraft  of  today  contain  an  INS  having  an  accuracy  of  approximately 
one  nautical  mile  (nm)  per  hour  when  operated  in  the  pure  inertial  unaided  mode.  A 
characteristic  of  this  class  of  INS  is  that  its  navigational  accuracy  is  determined  and  limited  by 
the  gyros  whose  technology  has  resulted  in  an  error  performance  of  approximately  0.0 1  degrees 
per  hour  drift  rate.  While  there  are  at  least  one  hundred  hardware  mechanizations  and  software 
computational  error  sources  that  contribute  to  the  navigation  errors,  these  are  generally 
controlled  such  that  the  gyro  drift  rate  errors  limit  INS  accuracy.  Figure  5.1  illustrates 


typical  INS  error  performance  for  a  local  level  platform  mechanization  having  a  north  gyro 
drift  rate  error  of  0.0 1  degrees  per  hour  ( 26:6 1 ).  These  characteristics  were  obtained  by 
simulating  a  set  of  ten  simultaneous  differential  equations.  The  error  performance  indicated  is 


Fig.  5. 1  Typical  INS  Navigation  Errors  (26:61 ) 


of  the  correct  order  of  magnitude  for  a  one  nm  INS.  Other  error  sources,  including  alignment, 
would  tend  to  modify  these  characteristics  but  would  also  result  in  platform  tilt  errors 
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exhibiting  the  84  minute  period  Schuler  oscillations.  An  exception  to  the  above  error  sources 

25 

would  be  th8  situation  of  violent  aircraft  maneuvers  coupled  with  acceleration  sensitive  gyro 
drift  rate  errors. 

j  Determining  the  feasibility  of  a  one  nm  per  hour  INS  for  performing  the  motion 

compensation  function  depends  to  a  large  extent  on  the  time  dependent  linear,  quadratic,  and 
higher  order  position  error  components.  The  "order  of  magnitude"  amplitude  of  these 
|  components  can  be  determined  by  expanding  the  longitude  error  of  Figure  5.1  into  a  Taylor 

series.  To  accomplish  this  expansion,  the  east  velocity  error,  Sv#,  is  approximated  by: 

Ve  =  cos(ust)  -  1  (ft/sec)  (5.1) 


c  • 

I 


where  t  is  time  in  seconds,  u>  is  2iT/Ts  and  Ts  is  the  Schuler  period  of  5040  seconds.  The  east 
position  ( longitude)  error  is  found  as  the  integral  of  the  velocity  error: 


P„  =  i(cosojt)  -  1 1  dt  =  - 1 


The  precision  of  the  position  information  AP#  rather  than  absolute  accuracy,  over  a  SAR 
coherent  integration  interval  AT ,  is  the  quantity  of  interest  here.  So,  j 

APe(AT)  =  Pe(t)  -  Pe(t  +  AT)  (ft)  (5.3) 

Substituting  equation  5.2  into  5.3  and  expanding  the  sin(t+  AT)  term  results  in 
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APe(AT)  =  AT  +  ^jsinwt-sinwtcoswAT-coswtsinwAT]  (5.4) 

Now  this  position  error  can  be  expanded  into  its  low  and  higher  order  time  components  by  the 
use  of  Tr/lor  series  expansions  and  specific  values  of  time: 

2  4  6 

coswAT  =  1  -  (WAT)  +  (waT)  -  (waT)  +  . . . .  (5.5) 

2!  4!  6! 

^ , . . . .  (5.6) 

to  arrive  at 


2  4  6 

(WAT)  +  (WAT) -(WAT)  +  . . . . 
2!  4!  6! 


(ft)  (5.7) 


In  a  similar  manner: 


3  5  7 

sinwAT  =  WAT  -(wAT)  +  (wAT)  -  (wAT)  + 


APe  =  — 

.  *  w 


3! 


GJt=  o,  21T  ,4ff,6lT, 


3  3  7 

(wAT)  -  (wAT)  ♦  (waT)  - 
3!  5!  7! 


7! 


Representative  values  of  these  time  dependent  errors  are  given  in  Table 
hour  INS  with  a  coherent  integration  time  of  3.5  seconds.  These  position  i 


(ft) 


(5.8) 


(5.9) 


(5.10) 


i.  1  for  a  one  nn»  per 
tot  components  are 


TABLE  5.1  REPRESENTATIVE  INS  POSITION  ERRORS 


Time  Dependency 

Aircraft  Motion 

Error  Magnitude  (ft) 

Linear 

Velocity 

-3.5 

Quadratic 

Acceleration 

00764 

Cubic 

Jerk 

.00001 

seen  to  drop  off  rapidly  with  increasing  order  and  it  is  this  very  characteristic  which  makes  an 
INS  a  feasible  motion  compensation  sensor  As  seen  previously,  SAR  image  quality  will  largely 
be  determined  by  the  quadratic  and  higher  order  components.  Table  5. 1  shows  that  these 
components  are  relatively  small  when  compared  to  a  wavelength  of .  1 06  feet  at  X-band 

Aided  INS 

While  it  has  been  shown  that  a  nominal  one  nm  per  hour  INS  has  an  inherent  capability  to 
f)  •  provide  the  desired  motion  sensing  for  many  S Aft  applications,  it  is  possible  that  there  are 

difficult  mapping  conditions  (long  range,  high  resolution,  aircraft  maneuvers,  etc.)  where 
higher  INS  performance  would  be  desirable  in  ordar  to  focus  the  SAR.  An  aided  inertial  Kalman 
filter  INS  mechanization  utilizing  velocity  and  position  updates  affords  a  means  of  greatly 
enhancing  INS  performance  with  regard  to  reducing  quadratic  phase  errors.  Table  5.2  provides 
an  assessment  of  state-of-the-art  sided  INS  performance  where  various  position  and/or 
velocity  measuring  sensors  have  been  Integrated  with  a  high  accuracy  0. 1  nm  per  hour  INS.  The 
system's  ability  to  accurately  measure  aircraft  acceleration,  velocity,  and  position  is  shown  and 
corresponds  to  a  rather  long  30  minuta  ground  alignment  anda  snort  20  mfnute  flight  time.  The 
radar  position  accuracy  given  reflects  INS  state  immediately  after  the  measurement.  A  top  level 
error  budget  showing  the  various  contributors  to  the  acceleration  measurement  errors  of  Table 
5.2  is  given  in  Table  5.3.  Data  is  given  for  straight  and  level  flight  and  a  6g  maneuver.  The 
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TABLE  5.2  PROJECT 

ED  AIDED  INS  PERFORMANCE  (16:1?) 

AIDED  INS  MECHANIZATION 

MEASUREMENT  ERROR 

Master-Slave  or 

Improved  Strapdovn 

Acceleration  (ug) 

Velocity 
(m  /sec) 

Position 

(meters) 

Level 

Flight 

6g 

Maneuver 

Unaided  System 

112 

715 

.3 

73 

Radar  Altimeter-DLMS  Aidinq 

112 

715 

.3 

75 

GPS  Aiding  (Type  1) 

50 

874 

06 

Velocity  Update  &  Star  Fix 

M2 

715 

.1 

73 

Unconv  God.  Align 

112 

500 

3 

75 

Radar  Position  Update 

112 

715 

.3 

12-23 

Current  State-of-the-Art 

Strapdown  (Laser  Gyros) 

Unaided  System 

13? 

3138 

i 

323 

Radar-DLMS  Aidino 

13? 

3158 

i 

75 

GPS  Aiding 

50 

874 

.06 

6 

Radar  Aiding  &  Star  Fix 

137 

715 

.6 

450 

Unconv  Gnd.  Align 

137 

3138 

1 

12-23 

accelerometer  null  shift  bias,  scale  factor,  nonorthogonality,  and  mechanical  misalignment 
errors  determine  the  residual  verticality  errors  after  ground  alignment  The  major  error  in 
the  gyrocompass  ground  alignment  Is  an  azimuth  misalignment  caused  by  the  east  component  of 
the  0.00 1  degree  per  hour  gyro  drift  rate.  It  is  important  to  point  out  that  gravity  anomalies 
affect  INS  accuracy  because  accelerometers  do  not  measure  acceleration  directly  They  measure 
specific  force,  which  is  the  vector  difference  between  acceleration  and  gravity.  Therefore,  it  is 
necessary  to  compute  the  gravity  vector,  and  to  add  It  to  the  vector  output  of  an  accelerometer 
triad  in  order  to  obtain  a  measure  of  acceleration.  The  computation  of  the  gravity  vector  entails 
the  usa  of  an  imperfect  mathematical  model  of  the  earth's  gravity  field.  The  errors  in  this  modal 
are  the  gravity  anomalies.  The  current  state-of-the-art  in  gravity  models  is  such  that  the 
anomalies  have  6.3  arc  sec  deflection  (errors  in  the  along- track  and  cross- track  direction), 
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TABLE  5  3  ACCELERATION  MEASUREMENT  ERROR  BUDGET  (16:18) 


ERROR 

CONTRIBUTION  TO 
MEASUREMENT  ERROR  (yg) 

1 .  Gr  avity  Anomalies 

tCT  Value 

6g 

Maneuvering 

a  RSS  Deflection 

9  arc  sec 

43 

43 

b  Magnitude 

68  ug 

68 

68 

2.  Accelerometer  Errors 

50  jig 

50 

50 

3.  RSS  Accelerometer  Platform  Drift 

4  Accelerometer  Platform  Ground 
Alignment 

0.001  deg/hr 

60 

10 

a.  RSS  Vertically 

(Accelerometer  Error) 

141  ug 

425 

71 

b  Azimuth  Misalignment 
(East  Gyro  Drift  Rate  During 
Gyrocompassing) 

0.001  deg/hr 

564 

0 

RSS 

715 

112 

i 

and  67.5  micro-g  magnitude  error  in  the  vertical  direction.  The  6.3  arc  sec  gravity  deflections 
are  equivalent  to  a  43  micro-g  acceleration  error  in  the  horizontal  plane.  Nothing  can  be  done 
*  at  present  about  this  kind  of  error  by  making  better  accelerometers  and  gyros.  • 


* 
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Vf.  SAR  notion  Compensation  Requirements 


It  has  been  shown  in  the  previous  sections  that  the  conceptual  compatibility  exists 
between  INS  performance  and  SAR  motion  compensation  requirements.  This  section  will  explore 
the  INS  performance  necessary  to  achieve  the  QPE  required  for  SAR  focusing  in  various  specific 
operating  conditions.  A  major  source  of  quadratic  phase  error  is  caused  by  an  uncompensated 
acceleration  along  the  line-of-sight  ( LOS)  to  the  map-point.  This  distance  error  is  the  product 
of  one  half  the  acceleration  error  and  time  squared  and  is  a  quadratic  function.  This  LOS 
acceleration  is  matte  up  of  two  components;  (1)  centripetal  acceleration  due  to  ideal  straight 
line  motion  of  the  radar  as  it  tlies  past  the  map  reference  point,  and  (2)  that  to  maneuvers 
or  disturbances  from  the  ideal  flight  path.  From  the  first  of  these  two  components,  an 
uncompensated  error  occurs  when  the  INS  has  a  velocity  error  component  normal  to  the  LOS  to 
the  map  reference  point.  An  uncompensated  error  from  the  second  component  results  when  tl« 
INS  has  LOS  acceleration  measurement  error. 

Although  these  QPE  ore  fundamental  to  INS  and  SAR  Integration,  there  are  other  sources  of 
QPE  which  can  arise.  It  can  be  shown  that  bistatic  SAR  incurs  a  rather  severe  QPE  due  to 
illuminator  and  receiver  position  errors  and  also  from  stable  oscillator  frequency  differences. 
Imperfect  radar  system  mechanizations,  of  which  there  are  many,  may  result  in  a  time  delay 
between  the  sensing  and  application  of  motion  data  during  aircraft  maneuvers  causing  a  very 
substantia!  QPE.  These  end  other  potential  sources  of  QPE  must  be  considered  in  a  system  level 
error  budget  for  any  given  SAR  design  but  will  not  be  covered  in  this  study. 

LflS  ACBPlfiCflliflQ  MBWMirammt  Error 

In  order  to  relate  the  required  INS  acceleration  measurement  accuracy  to  SAR  azimuth 
resolution,  aircraft  motion,  and  scenario  parameters,  an  equation  is  derived  for  the  LOS 
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acceleration  measurement  accuracy.  The  INS  requirement  which  results  is  then  compared  to 
current  INS  performance  capabilities  to  determine  the  nominal  SAft  operating  boundaries 
imposed  by  this  particular  motion  compensation  error.  In  order  to  get  the  equation  for  the  LOS 
acceleration  measurement  accuracy  the  quadratic  term  of  the  Taylor  scries  expansion  in 
Equation  42  is  substituted  into  the  expression  for  phase  error  as  given  by  Equation  4.3. 

A4/q(t)=  —  jsr  (radians)  (6.1) 

A4^(T/2)=  ^A8losT"  (radians;  (6.2) 

where: 

AaL0S§i^l)  md  t  =  I  (6.3) 

If  Equation  2.22  is  rearranged  to  obtain  the  coherent  integration  time  and  a  permissible 
quadratic  phase  error  of  90  degrees  Is  used: 

2 

A%(T/2>  =  AO L0S  zl  (radians)  (6.4) 

8dzVzsinz0  2 

4d2V2s?n20 

aaios  = - 5 -  (ft/sec)  (6.5) 

AR2 

It  must  be  remembered  that  this  expression  for  the  LOS  acceleration  measurement  accuracy  is 
optimistic  in  that  it  requires  the  entire  error  budget  and  makes  no  allowance  for  a  centripetal 
acceleration  or  other  OPE  components.  Required  acceleration  measurement  accuracy  vs.  squint 
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angle  is  tabulated  in  Table  6. 1  to  show  the  increasing  accuracy  needed  as  the  squint  angle 
decreases.  Also,  in  Table  6.2,  the  acceleration  measurement  accuracy  vs.  tactical  aircraft  speed 
is  tabulated  showing  the  increasing  accuracy  needed  as  aircraft  speed  decreases.  A  comparison  of 
the  motion  compensation  requirements  of  Tables  6. 1  and  6.2  and  the  acceleration  measurement 
accuracy  performance  of  Table  5.2  indicates  that  a  state-of-the-art,  0.1  nm  per  hour,  aided  INS 


TABLE  6.1  LOS  ACCELERATION  ACCURACY  REQUIREMENTS 
VS.  SQUINT  ANGLE 


^aLOS<US> 

Velocity  (ft /see) 

760  (450  knots) 

912  (540  knots) 

Azimuth 

Squint  Angie  (  °) 

Squint  Angle  (  °) 

Res  (m) 

20  30  40  50  90 

20  30  40  50  90 

1 

88  1SS  311  441  752 

127  271  447  636  1083 

3 

792  1692  2797  3972  6769 

1140  2437  4027  5720  9747 

10 

8798  18803  31075  44136  75211 

12669  77076  44749  63556  108304 

R=30  km=9.84252  x  1 0**  ft  A*. 106  ft 

TABLE  6  2  LOS  ACCELERATION  ACCURACY  REQUIREMENTS 
VS.  AIRCRAFT  SPEED 


AaL0S(»9> 

7€n  (450  knots) 

Velocity  (ft/sec) 

810  (480  knots) 

912  (540  knots) 

Azimuth 
Res  (m) 

Squint  Angle  (  °) 

20  30  90 

Squint  Angle  (  °) 

20  30  90 

Squint  Angle  ( °) 

20  .  30  90 

1 

88  188  752 

100  214  854 

127  271  1083 

3 

792  1692  6769 

899  1922  7689 

1140  2437  9747 

10 

8798  18803  75211 

9994  21358  85433 

12669  27076  108304 

R=30  km=9 .84252' x  10*  ft  A  =  .  1 06  ft 


could  adequately  motion  compensate  a  SAR  with  resolution  as  fine  as  3  meters  under  most 

! .  .  . 
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operating  conditions.  Only  during  high  g  air  craft  maneuvering  and  smaU  squint  angie  operation 
does  the  SAR  requirement  become  exceeded  by  the  INS  capability.  Achieving  even  finer 
resolutions  than  three  meters  would  be  possible  if  imaging  was  done  at  ranges  shorter  than  30 
km. 

LOS  Centripetal  Arrfiiffration  Error 

If  the  geometry  in  Figure  6. 1  is  considered,  the  velocity  accuracy  required  to  limit  the 
LOS  centripetal  acceleration  error  to  a  level  consistent  with  SAR  OPE  can  be  determined.  The 
instantaneous  aircraft  range  to  the  map  reference  point  can  be  expressed  as: 

R(0  =  [rJ+  (Vt)2]2  (ft)  (6.6) 

where  R0  corresponds  to  the  time,  t=0 ,  at  which  the  aircraft  comes  closest  to  the  map  reference 
point 


FIG.  6.1  Aircraft  Closure  Geometry 
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Velocity  accuracy  requirements  are  shown  in  Table  6.3  and  should  be  compared  to  the  projected 
aided  INS  performance  capabilities  shown  in  Table  5.2  of  roughly  one  ft/sec.  Again  it  is  shown 
that  state-of-the-art  INS  are  capable  of  providing  motion  compensation  for  SAR  resolutions  as 
fine  as  three  meters  for  a  wide  range  of  scenarios. 
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VII.  Autofocus  Techniques 


It  has  been  shown  that  a  high  degree  of  compatibility  exists  between  SAR  motion 
compensation  requirements  and  INS  performance  capabilities  but  high  resolution  imaging  at  long 
range,  small  squint  angles,  low  aircraft  velocities,  and  high  g  maneuvers  may  impose  QPE 
requirements  which  cannot  be  satisfied  by  aided  INS  cap^iilities.  In  such  situations,  it  may  be 
necessary  to  perform  a  post-collection  motion  compensation  correction  prior  to  the  Fast 
Fourier  Transform  (FFT).  It  may  be  noted  that  normal  motion  compensation  takes  place 
"on- the- fly"  as  the  return  d3ta  is  received  by  the  radar  or  as  a  closely  associated  "one-pass" 
process.  Autofocusing  is  a  technique  which  extracts  information  from  the  partially  processed 
image  in  order  to  estimate  the  residual  error  phase  which  is  then  taken  out  of  the  date  prior  to 
the  full  processing.  Figure  7. 1  illustrates  the  interface  between  the  autofocusing  process  and 
the  azimuth  processing.  As  can  be  seen,  the  two  azimuth  processing  steps  are  an  essential 
prerequisite  for  many  types  of  autofocus  techniques. 


Fig.  7.1  Autofocus  Interface  (3:  484) 


in  the  last  few  years,  several  approaches  have  been  proposed  to  estimate  the  amount  of  QPE 
present  in  a  radar  signal.  Some  of  the  earlier  methods  require  certain  map  structures  to 


perform  successfully,  ie,  they  require  the  existence  in  the  map  of  strong  isolated  point 
scat  ter  ers.  It  would  he  desirable  to  obtain  an  autofocus  algorithm  that  performs  on  any  scene 
structure.  An  algorithm  to  accomplish  this  is  called  the  map  drift  method  of  autofocus. 

hap  Drift  Method 

The  utilization  of  map  drift  to  measure  SAR  misfocus  dates  back  to  the  observation  by 
operators  of  early  SAR  systems  that  focus  errors  cause  azimuth  drift  of  successive  maps  of  the 
same  terrain.  This  drift  is  caused  by  quadratic  phase  variation  in  the  raw  data,  which,  when 
integrated  over  the  time  difference  between  two  maps  of  the  same  scene,  results  in  a  linear 
(with  time)  phase  difference  between  the  data  of  the  two  maps.  The  apparent  Doppler  shift 
caused  an  azimuth  offset  between  the  two  maps.  This  source  of  map  drift  can  be  illustrated  by 
Figure  7.2.  A  quadratic  phase  variation  .assumed  zero  at  the  center  of  the  array,  is  illustrated 


Fig.  7  2  Subarray  Linear  Phases  with  Opposite  Slopes 
Caused  by  Quadratic  Phase  on  Full  Array  (3:  486) 


If  the  array  is  divided  into  two  subarrs/s,  the  quadratic  phase  an  each  side  can  be  considered  to 


consist  of  quadratic  variation  whose  maximum  value  is  one  fourth  that  of  the  full  array,  in 
addition  to  a  linear  phase  term.  The  difference  in  the  linear  phase  terms  (constant  frequency 
difference)  causes  the  azimuth  map  offset  The  drift  between  subarray  maps  can  be  measured 
several  different  ways  and  each  is  a  candidate  for  implementation  as  a  quadratic  phase  estimator. 
The  most  obvious  method  is  to  cross  correlate  the  subarray  maps  and  choose  the  lag  or  offset  that 
maximizes  the  correlation  integral.  This  approach  is  outlined  in  Figure  7.3 .  The  input  is  a 


Fig.  7.3  Block  Diagram  of  Map  Drift  Aitofocus  Algorithm  ( 1 6:25) 


sequence  of  data  that  has  been  focused  approximately  by  using  the  best  knowledge  of  the  two 
radar's  velocity,  acceleration,  and  antenna  look  angles.  The  autofocus  algorithm  operates  on  a 
subarray  of  data,  of  width  leas  than  1  /2  of  the  full-reaolution  array,  and  forms  a  map  from  the 
magnitudes  of  the  Fourier  components  after  Fourier  transformation.  The  process  is  repeated 
later  with  some  nonoverlapping  portion  of  the  phase  history  of  the  same  ground  scene  at  the 
same  range.  Tha  two  maps  are  cross-correlated,  and  the  resulting  cross-cnrrelation  function  is 
averaged  over  range  to  smooth  out  scintillation  and  other  noise  disturbances.  The  estimated  drift 
in  pixels  is  then  taken  to  be  the  offset  at  which  the  averaged  cross  correlation  function  peaks. 
The  accuracy  of  the  offset  can  be  estimated  to  within  less  than  one  pixel  by  using  an  interpolating 
function  to  find  the  peak  of  the  sampled  cross  correlation  function.  When  the  subarrays  are 
spaced  one- half  of  the  full  resolution  array  length  N,  the  offset  in  number  of  pixels  on  the 
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subarray  map  multiplied  by  it  is  the  estimated  residual  QPE  in  radians  measured  center  to  end  of 
the  full  array.  Cnee  the  QPE  has  been  determined,  the  full  array  map  data  is  refocused  and  then 
the  map  is  Fast  Fourier  Transformed  and  displayed. 

SAP  Performance  Envelopes 

As  illustrated,  the  worst  case  motion  compensation  conditions  indicate  difficult 
requirements  on  velocity  and  acceleration  motion  sensing  equipment  If  sensing  and  compensation 
were  tn  be  performed  in  one  of  the  traditional  ways.  The  projected  state-of-the  art  capability  in 
motion  sensing  techniques  is  on  the  order  of  50  to  100  pg  acceleration  and  0.2  to  0.5  m/sec 
velocity  when  considered  in  terr.  s  of  operation  in  a  highly  maneuverable  aircraft. 
Consequently,  a  reliance  on  traditional  motion  compensation  techniques  as  the  only  approach  to 
the  motion  compensation  problem  would  result  in  degraded  performance  at  the  extremes  of  the 
operational  envelope.  Autofocus  technology  has  progressed  significantly  in  the  past  few  years  to 
»  the  point  it  can  be  exploited  to  complement  traditional  motion  compensation  approaches  to 

significantly  extend  the  conditions  over  which  desired  performance  can  be  achieved. 
Assessments  of  SAR  have  been  made  on  various  programs  with  the  conclusion  that  autofocus 
techniques  are  well  suited  to  accommodating  large  QPE.  The.  upper  bound  on  resolution 
improvement  which  can  be  provided  by  autofocus  is  a  complex  function  of  clutter-to- noise 
ratio,  scene  content,  scene  size,  and  higher  order  phase  errors  end  Is  just  now  becoming  a  topic 
of  both  theoretical  and  experimental  investigation.  QPE  in  the  region  of  24ir,  32ir,  and  40tt, 
have  routinely  been  accurately  estimated  and  compensated  in  recent  programs  utilizing  high 
clutter-to- noise  ratio  data  with  artificially  introduced  QPE.  ThSeQPE  represent  nominal  values 
and  the  maximum  achievable  QPE  compensation  has  not  yet  been  determined 

To  conclude  this  study,  Figure  7.4  illustrates  the  significance  of  various  levels  of  motion 
compensation  with  and  without  autofocus.  The  graphs  in  this  fllustrotion  show  the  rsngs-squint 
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angle  regions  where  focus  can  be  maintained  based  on  the  indicated  velocity  and  acceleration 
compensation  capabilites.  The  lower  set  of  graphs  indicate  the  performance  envelope  which  can 
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Fig.  7.4  Performance  Envelope  Sensitivity  to  Various 
Levels  of  Motion  Compensation  Capabilty  ( 1 6:27) 

be  achieved  without  autofocus  while  the  above  set  shows  the  capability  that  can  be  achieved  with 
a  1  Ox  improvement  via  autofocus.  The  1 0x  value  is  a  baseline  for  comparison  and  was  selected 
to  reflect  a  system  mechanization  where  the  inertial  motion  compensation  is  designed  to 
accomodate  a  10m  moderate  resolution  capability  in  the  first  FFT  process  with  a  tm  fine 
resolution  capability  being  provided  by  a  second  FFT  process.  The  I  Ox  relates  to  the  fact  that 
the  ration  of  moderate  to  fine  resolution  is  10.  This  improvement  value  assumes  a  OPE 
reduction  of  5 On  by  the  autofocus  algorithm.  The  outer  squint  angle  boundaries  are  determined 
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by  somewhat  arbitrary  radar  waveform  ambiguity  constraints  associated  with  second  time 
around  range  echoes  and  the  higher  pulse  repetition  frequencies  necessary  when  imaging  at  the 
larger  squint  angles.  Figure  7.4  indicates  that  the  operational  envelope  is  significantly 
compromised  even  with  a  good  ( lOOug  and  1  m/sec)  motion  compensation  capability  without 
autofocus.  On  the  other  hand,  with  a  lOx  autofocus  improvement,  a  very  good  performance 
envelope  is  achievable  with  what  could  be  termed  as  a  poor  ( I  mg  and  3  m/sec)  motion  sensing 
capability.  Figure  7.5  shows  the  ability  of  autofocus  augmentation  to  maintain  the  range-squint 
angle  performance  envelope  in  the  presence  of  INS  acceleration  accuracy  degradations  during 
aircraft  high  g  maneuvers. 


Motion  Comp  Constraint— INS  Aufofocus 
(100  ^g/0.3  m/sec  *  10*  Autofocus) 

1g  3g  6g 


F ig.  7.5  Aircraft  Maneuver  Effects  on  EAR  Range-Squint  Angle  Focus  Reqions 


VIII.  Conclusions  and  Recommendations 


This  study  has  demonstrated  the  feasibility  of  using  state-of-the-art  strapdown  inertial 
Navigation  Systems  to  satisfy  the  motion  compensation  requirements  for  Synthetic  Aperture 
Radar  onboard  maneuvering  aircraft.  This  was  accomplished  by  deriving  the  equations  which 
related  Inertial  Navigation  System  capabilities  to  Synthetic  Aperture  Radar  motion  compensation 
requirements.  Both  tuned  and  ring  laser  gyro  systems  were  investigated  to  determine  their 
applicability  for  the  task.  In  particular,  the  autofocus  technique  was  discussed  for  motion 
compensation.  Autofocus  is  one  means  of  performing  the  motion  compensation  function  aboard 
aircraft  that  ore  maneuvering  os  apposed  to  benign  straight  and  level  flight.  Autofocus  was 
shown  to  improve  the  range  and  squint-angle  boundaries  significantly.  For  a  moderate  lOx 
improvement  via  autofocus,  it  was  shown  that  the  performance  envelope  for  a  good  ( 1 0Opg  and  1 
m/sac)  motion  compensation  capability  was  extended  considerably.  Far  a  poor  ( tmg  and 
3m/sec)  motion  compensation  capability,  it  was  possible  with  outofocus  to  have  a  very  useful 
performance  envelope 

The  upper  bound  on  resolution  improvement  which  can  be  provided  by  autofocus  is  a 
complex  function  of  clutter-to- noise,  scene  content,  scene  size,  and  higher  order  phase  errors 
and  should  be  investigated  further  to  determine  the  areas  in  which  study  should  be  accelerated.  A 
further  investigation  into  this  complex  function  should  be  undertaken  to  determine  which  areas 
provide  significant  improvement  in  resolution.  Also,  the  outer  squint  angle  boundaries  have 
bean  determined  using  arbitrary  radar  waveform  ambiguity  constr  aints  and  these  need  to  be 
investigated  further  to  ascertain  their  appropriateness. 
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